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1.0  IDENTIFICATION 

1.1  Compute  shadow  times  (SHADOW).  IA5  MOD  AA. 

1.2  By  C.  M.  Seacat,  15  May  1963 

1.3  System  Development  Corporation,  Sunnyvale,  California 
2.0  PURPOSE 

SHADOW  is  a  l60l  COP -controlled  program  to  compute  the  times  at  which  a 
satellite  enters  and  leaves  the  shadow  of  the  earth.  Results  are  printed  on¬ 
line  and  off-line. 

3-0  USAGE 


3  •  1  Calling  Sequence 

The  program  is  normally  called  by  a  function  card  with  the  format: 


"SHADOW  A  B  C  D  E  F.F  G  H  I  J  K.K 


A  Vehicle  number 
B  Month 

C  Day 

D  Hour 

E  Minute 

F.F  Second  (Floating  Point) 

G  Month 

H  Day 

I  Hour 

J  Minute 

K.K  Second  (Floating  Point) 


Start  Time 


Stop  Time 
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An  option  exists  to  specify  start  and  stop  times  by  rev  number. 
The  function  card  format  is: 


♦SHADOW  A  B.  C. 


A  Vehicle  Number 

B  Starting  Rev  (decimal  point  mandatory) 

C  Ending  Rev  (decimal  point  mandatory) 

A  fractional  rev  number  may  be  input  without  harming  the  program, 
but  the  fractional  part  will  be  ignored.  For  consistency  between 
the  function  card  and  printed  output,  it  is  recommended  that  zero, 
or  blank,  follow  the  decimal  point  in  the  rev  number.  Calculations 
are  carried  out  to,  but  not  including, the  terminal  rev.  In  this  option 
starting  rev  must  be  1.  or  greater. 

Reset  tapes  containing  orbital  conditions  must  be  mounted  on  tape 
drives  2,  9  and  10. 

There  are  no  unexplained  halts  in  the  program  with  one  exception. 

Early  in  the  program  the  subroutine  SETOP  is  called.  If  reset  tape 
and  reference  pool  parameters  are  not  correct  SETOP  will  halt  due 
to  an  erroneous  return  Jump.  The  trouble  1b  likely  to  be  in  the 
timing  area,  for  example  a  launch  time  earlier  than  MONTHB  in  the 
reference  pool. 


3-2  OUTPUTS 

Every  time  the  satellite  leaves  the  shadow  of  the  earth,  the  program 
prints  on  line  the  date,  time,  rev  number,  geodetic  latitude,  and 
longitude  at  which  it  entered  and  left.  When  the  stop  time  is 
reached,  the  program  prints  out  the  percentage  of  time  in  shadow, 
i.  e.,  the  sum  of  all  the  shaded  periods  divided  by  the  time  range 
on  the  function  card. 

If  the  satellite  does  not  enter  shadow  for  a  particular  rev  this 
Information  is  printed. 

If  the  orbit  decays  during  the  time  span  of  the  function  card  this 
Information  is  printed. 

All  output  is  on-line  and  duplicated  off-line  (tape  3).  Figure  1 
shows  sample  printouts. 
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FIGURE  1. 
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To  start  from  scratch,  the  following  function  sequence  Is  recommended: 
*IRT  ... 

*WNRT  . . .  (Bach  card  completed  appropriately) 

♦HUIUI  ... 

♦SHADOW  ... 


4.0  METHOD 

4.1  Mathematical  Formulation 

The  following  assumptions  are  made: 

a.  the  shadow  of  the  earth  is  cylindrical  and  non-fuzzy . 

b.  the  positions  of  the  satellite  and  the  sun  are  known 
for  a  given  time,  t,  in  inertial  coordinates. 

The  coordinates  of  the  sun  at  time  t  are  (x,  y,  z).  The  axis  of 
the  cylinder  formed  by  the  shadow  of  the  earth  passes  through  the 
point  (-X,  -y,  -z). 

The  coordinates  of  the  satellite  at  time  t  are  ( 5, T)  ,i  ). 

The  equation  p2-  52+Tp  +£2  defines  an  earth-centered  sphere  of 
radius  p,  with  the  satellite  on  the  surface. 

The  axis  of  the  cylindrical  shadow  pierces  this  sphere  at 
(x1 ,  y' ,  z' ),  where: 

*’  '  --P 

y'  ■  -  _£l 
R 

z '  *  -  Pz 

R 


V2  2  2 

x  +  y  +  z  ,  distance  of  sun  at  time  t 
The  geometry  is  illustrated  in  Figure  2. 

The  distance,  d,  from  (x1,  y' ,  z' )  to  the  circular  intersection  of 
the  cylinder  and  the  sphere  is  given  by: 


d2  -  2  P  {  P  -  \|p  2  -  r2 


)  where  r  ■  earth  radius. 
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The  distance,  D,  from  (x',  y' ,  z1 )  to  the  satellite  is  given  by 

D2  «  (x*  -if  +  (y'  -Tl)2  +  (z'  -C)2 

If  D2>d2,  the  satellite  is  in  sunlight. 

2  2 

If  D  *d  >  the  satellite  is  in  shadow. 

Please  note  that  the  analysis  above  does  not  assume  a  circular 
or  near-circular  orbit.  A  new  ^  is  computed  for  each  t. 

4.2  Coordinate  Determination 

The  program  computes  the  inertial  coordinates  of  the  satellite  at 
time  t  by  the  subroutines  TTE  and  PTR. 

The  coordinates  of  the  sun  are  computed  by  the  program  using  an 
approximation  technique  based  on  the  near-sinusoidal  variation  in 
sun  coordinates  as  a  function  of  time.  The  crucial  parameters  are 
the  equlnoces  (at  which  time  y  and  z  are  zero),  the  times  at  which 
y  and  z  are  maximum  and  minimum,  the  times  at  which  x  is  ma-irt  mum 
and  minimum,  and  the  respective  maximum  and  minimum  values.  Figure  3 
is  a  plot  of  sun  x,  y  and  z  against  time. 

The  trick  is  to  fit  x,  y,  and  z  with  sine/cosine  curves  as  nearly 
as  possible.  Bach  curve  in  figure  1.  is  broken  into  three  parts: 
pre -vernal  equinox,  between  equlnoces,  and  post-autumnal  equinox. 

(in  the  case  of  x,  the  schisms  are  not  exactly  at  the  equlnoces, 
but  at  the  maximum  and  minimum  values,  which  are  close  to  the 
equlnoces).  A  slightly  different  period  is  used  for  each  case,  the 
choice  depending  on  where  in  the  calendar  year  t  falls.  The  amplitude 
of  y  and  z  are  slightly  different  outside  the  equlnoces  than  between 
them,  for  which  the  program  condensates .  Finally,  x  is  adjusted 
downward  by  a  constant  amount.  Examination  of  the  table  below  and 
the  pertinent  sections  of  the  flow  charts  will  clarify  the  details 
of  this  technique. 

The  program  will  be  valid  as  is  through  1972.  At  that  time  the 
parameters  tabulated  below  should  be  re -computed  and  the  tagged 


locations 

updated. 

TAG 

UNITS 

SCALING 

DESCRIPTION 

XBASE 

Hours 

Floating  Point 

Time  at  which  x  is  maximum 

XBASE 

Hours 

Floating  Point 

Vernal  equinox  (y,  z  ■  0) 

XPER 

Hours 

Floating  Point 

Difference  between  times  at 
which  x  is  maximum  and  minimum. 

XPER 

Hours 

Floating  Point 

Difference  between  equlnoces. 
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FIGURE  3.  SUN  COORDINATES  VERSUS  TIME 
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TAG 

XPERN 

UNITS 

Hours 

SCALING 

Floating  Point 

XPERN1 

Hours 

Floating  Point 

YPERN 

Hours 

Floating  Point 

YPERN1 

Hours 

Floating  Point 

XADJ 

None 

Floating  Point 

YNORMP 

None 

Floating  Point 

YNOHMN 

None 

Floating  Point 

ZNOHMP 

None 

Floating  Point 

ZNOKMN 

None 

Floating  Point 

THISYR 

Integer 

Right  Just 

XBASES 

Hours 

Floating  Point 

XBASES+1 

Hours 

Floating  Point 

DESCRIPTION 

Twice  the  Interval  between 
XBASE  and  the  time  In  the 
preceding  year  at  which  x  Is 
equal  to  XADJ,  defined  below. 
Ibis  time  will  be  close  to  the 
last  year's  winter  solstice. 

Twice  the  interval  between  the 
time  at  which  x  is  minimum  and 
the  time  this  year  that  x  equals 
XADJ,  defined  below.  This  time 
will  be  near  this  year's  winter 
solstice. 

Twice  the  interval  between  the 
time  at  which  y  is  minimum  the 
preceding  year  (near  winter 
solstice)  and  this  year's 
vernal  equinox. 

Twice  the  interval  between 
this  year's  autumnal  equinox 
and  the  time  at  which  y  is 
minimum  this  year  (near  winter 
solstice). 

Minus  one  half  the  sum  of  the 
maximum  and  minimum  values  of  x. 

The  maximum  value  of  y. 

The  absolute  value  of  the 
minimum  value  of  y. 

Hie  maximum  value  of  z . 

The  absolute  value  of  the 
minimum  value  of  z. 

Hie  year  for  which  these 
calculations  are  being  made. 

XBASE  for  this  year 

XBASE  +6.15  hours  for  next 
year. 


4 


{ 
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TAG 

UNITS 

SCALING 

DESCRIPTION 

XBASES+2 

Hours 

Floating  Point 

last  entry  plus  6.15  hours 

• 

for  next  year 

• 

• 

etc 

• 

• 

• 

• 

« 

etc 

• 

• 

XBASES+9 

• 

Hours 

Floating  Point 

• 

Note:  Subtract  24  from 
entries  corresponding  to 
year  after  leap  year. 

YBASES 

Hours 

Floating  Point 

YBASE  for  this  year 

YBASES+1 

• 

Hours 

• 

Floating  Point 

• 

YBASE  +  6.15  hours 

• 

etc 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

YBASES+9 

• 

Hours 

• 

Floating  Point 

• 

(Same  as  above,  for  Y) 

The  inertial  coordinates  of  the  sun  for  a  given  year  are  tabulated 
in  the  American  Ephemerls  and  Nautical  Almanac.  This  volume  should 
be  used  to  compute  the  parameters  required.  YNORMP,  YNOHMN,  ZNORMP, 
and  ZNORMN  are  transcribed  directly  from  the  almanac.  Appropriate 
conversion  is  done  by  the  program. 

4.3  Detailed  Program  Processes 

When  first  entered,  the  program  searches  for  the  boundaries  of  the 
shadow.  This  is  done  by  incrementing  TAU  by  one-eighth  of  a  period, 
computing  orbital  elements,  and  applying  the  equations  of  section  4.1. 
When  the  satellite  is  found  to  be  in  shadow,  the  search  interval  is 
reduced  to  one  minute,  and  the  orbit  advanced  from  the  last  knovn  time 
it  was  in  sunlight  (one -eighth  period  ago)  in  one  minute  increments 
until  it  again  enters  shadow.  The  corresponding  time  and  position 
axe  saved.  The  orbit  is  advanced  in  one-minute  intervals  until  the 
satellite  enters  sunshine.  One  line  of  data  is  printed. 

If  the  satellite  was  in  shadow  at  the  starting  time  of  the  function 
card,  the  "search"  logic  Just  described  is  repeated,  starting  at  the 
time  the  satellite  entered  sunshine.  Thus  after  one  or  two  revs  the 
shadow  boundaries  have  been  determined  to  within  one  minute.  Accordingly, 
the  first  and  perhaps  the  second  line  of  printout  show  one-minute 
(  precision. 
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After  the  boundaries  are  determined,  the  orbit  is  updated  to  last 
time  of  entry  into  shadow  plus  the  orbital  period.  IAU  is  advanced 
or  regressed  in  one-minute  increments  until  the  start  of  shadow  is 
bracketed,  then  advanced  from  sunlight  to  shadow  in  one-second 
Intervals.  The  time  of  entry  is  again  saved.  This  time  is  increased 
by  the  length  of  time  in  shadow  last  rev,  end  of  shadow  bracketed 
to  within  a  minute,  and  finally  resolved  to  one  second.  One  line 
of  data  is  printed  out.  This  and  subsequent  printouts  show  one 
second  precision.  This  is  repeated  until  IAU  exceeds  stop  time  from 
the  function  card.  This  logic  is  based  on  the  fact  that  due  to 
the  slow  change  in  sun  position,  successive  entering  times  are 
close  to  a  period  apart.  However,  sun  position  i£  recomputed  for 
each  TAU. 

If  the  shadow  interval  is  less  than  one  minute  for  a  particular  rev, 
this  is  printed  out  for  each  rev  that  it  occurs.  After  the  first 
occurrence  of  total  sunshine  the  search  logic  is  repeated  until  shadow 
Intervals  exceed  one  minute.  No  harm  is  done  if  the  satellite  starts 
in  total  sunshine. 

Detailed  flow  charts  may  be  found  in  section  10.0. 

5.0  RESTRICTIONS 
5.1  Accuracy 

It  was  noted  in  the  preceding  section  that  once  the  program  got  off 
the  ground  results  were  printed  to  the  nearest  second.  Bear  in 
mind  that  this  is  merely  the  precision  to  which  the  program  computes 
and  should  not  be  construed  as  a  claim  of  one  second  accuracy. 

Regarding  accuracy,  the  two  major  error  contributors  are  the 
approximation  of  sun  coordinates  by  sine  curves,  and  the  assumption 
that  the  earth's  shadow  is  a  cylinder  without  penumbra.  (Erroneous 
orbited  prediction  say  be  discounted). 

The  maximum  relative  error  due  to  sun  position  approximation  is 

+  .001. 

The  earth  umbra  is  of  course  conical.  The  ratio  of  height  to 
base  radius  is  about  410  to  1,  suggesting  a  relative  error  increase 
of  about  .002  for  each  earth  radius  the  satellite  is  distant. 
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The  effect  of  the  penumbra  will  not  be  considered  except  to  point 
out  that  it  would  tend  to  make  some  arbitrarily  defined  "mean 
shadow"  (e.g. ,  half-way  into  the  penumbra)  more  nearly  cylindrical. 
Refraction,  however,  works  against  this,  having  the  effect  of  shorten¬ 
ing  the  earth  shadow  core. 

These  error  limits  are  well  within  tolerances  of  expected  users 
of  the  program. 

5.2  Hardware  Requirements 

One  1604  computer  plus  normal  peripheral  equipment. 

5.3  Subroutines 

The  program  uses  the  following  subroutines. 


Code  Name 

Code  Name 

RESET 

FIX 

PTR 

RALONG 

LTT 

UTMAC 

TTE 

MACGUT 

SEHJP 

SORT 

SUBERR 

ATAN1 

RADDEG 

COS 

OUTPUT 

SIN 

FLOAT 

TAN 

The  program  uses  the  Breakwell  closed  form  equations,  ibis  introduces 
some  (not  well  defined)  restrictions  as  to  admissible  orbited 
eccentricity  and  nagnitude  of  semi -major  axis.  Pending  further 
study,  all  that  can  be  said  is  try  it  and  see.  Difficulties,  if 
any,  will  occur  in  TTE. 
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5.4  Reference  Pool 

The  program  refers  to  the  following  items; 


Symbol 

Set 

Used 

BIATE 

X 

BWPER 

X 

CLOCK 

X 

X 

MY 

X 

X 

MYTAB 

X 

DBG 

X 

X 

HRS 

X 

X 

IDBCL 

X 

IAT 

X 

X 

MIN 

X 

X 

MONTH 

X 

X 

RA 

X 

X 

RAD 

X 

X 

RE 

X 

REV 

X 

X 

SEC 

X 

X 

TAU 

X 

X 

X 

X 

X 

Y 

X 

X 

YEAKB 

X 

Z 

X 

X 

6.0  TIMING 

The  program  requires  at  most  three  seconds  per  rev,  and  about  one  second 
per  rev  average,  including  output  time. 
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7.0  STORAGE  REQUIREMENTS 

The  program  SHADOW  occupies  470  cells,  plus  42  spares  for  512  (lOOOg). 
Subroutine  requirements  can  be  ascertained  from  the  following  storage 
analysis: 


STORAGE  ANALYSIS 

(*  =  RELOCATABLE  N  -  NEWLY  DEFINED  R  -  REDEFINITION) 


Name 

Address 

Name 

Address 

RESET 

*  11000 

UTMAC 

*  21232 

TIE 

*  21541 

FLQATBIN 

*  22726 

RADDEG 

*  22544 

ATAN1 

*  22460 

FIX 

*  21203 

SIN 

*  22252 

MACTIM 

*  22630 

LTT 

*  21252 

UNPACK 

*  22714 

SUBERR 

*  20520 

ATAN2 

*  22744 

FLOAT 

*  22154 

COS 

*  22167 

RALONG 

*  22573 

SHADOW 

*  10000 

MACGUT 

*  21443 

PTR 

*  22067 

SORT 

*  22336 

SETUP 

*  17123 

AS  IN 

*  23057 

OUTPUT 

*  17176 

TAN 

*  22375 

REDLO  *  22604 

8.0  VALIDATION  AND  TEST  CASES 

In  order  to  verify  program  outputs  the  following  tests  were  performed. 

Function  cards  and  pertinent  computer  output  sure  shown  with  each  test. 

8.1  Test  1.  The  injection  vector  was  chosen  so  that  the  orbit  be  as 
nearly  circular  as  possible,  be  perfectly  polar,  contain 
the  sun  in  its  plane  at  the  summer  solstice,  and  have  a 
zero  drag  factor.  Under  these  conditions  the  shadow 
interval  should  be  maximum  at  the  summer  and  winter 
solstice.  Furthermore,  it  is  possible  to  examine  the 
cross  section  of  the  earth  and  satellite  orbit  cut  by  the 
orbital  plane  to  check  the  entering  and  leaving  latitudes 
output  by  the  program.  Figures  4  and  5  show  the  results . 
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FIGURE  5- 


'CROSS  SECTION,  ORBITAL  PLANE  CONTAINING  SUN  AT  WINTER  SOISTICE 
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Function  Cards  and  Results 


*  WNRT  1136  0  0 

R  37  DEC  0 

R  46  DEC  0 

R  60  DEC  0 

END 

*  TIME  1136  1  1  6  20  19  42  20.95 

*  TIME  1136  1  2  6  20  19  44-  55.58 

*  TIME  1136  1  3  6  20  19  45  32.28 

*  TIME  1136  1  4  6  20  19  49  36. 

*  injiun  1136  0  -25525.  -270.  +22.7  +117*4  -692250.  0.0  0 

*  CLR 

*  CARDS 

DEFINE  SHADOW  3000 
(BINARY  DECK) 

*  SHADOW  1136  6  20  20  0  0.0  6  24  0  0  0.0 


VEH  1136  RUN  NO.  TST11MAY 
SHADOW  INTERVALS  TABULATED  BELOW 


ENTER 


LEAVE 


MO  DY  HR  MN 


6  21  5  7 
6  21  6  35 
6  21  8  4 
6  21  9  33 
6  21  11  1 


6  23  5  52 
6  23  7  21 
6  23  8  50 
6  23  10  18 
6  23  11  47 


SEC 

SYSTEM 

REV 

HR  MN  SEC  SYSTEM 

• 

REV 

13. 

18433- 

6.7 

• 

5  44  7-  20647. 

7-1 

52. 

23752 

7-7 

7  12  46.  25966. 

8.1 

32. 

29072. 

8.7 

8  4l  26.  31286. 

9-1 

11. 

34391- 

9-7 

10  10  5-  36605. 

10.1 

51. 

39711- 

10.7 

11  38  45.  41925. 

• 

11.1 

56. 

21176. 

39-7 

• 

6  29  50.  23390. 

40.1 

35* 

26495 

40.7 

7  58  28.  28708. 

41.1 

15 

31815 

41.7 

9  27  8.  34028. 

42.1 

54. 

37134. 

42.7 

10  55  47-  39347- 

43.I 

34. 

42454. 

43-7 

12  24  27.  44667. 

44.1 

DURATION 


2214. 

2214. 

2214. 

2214. 

2214. 


2214. 

2213. 

2213. 

2213. 

2213. 
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*  SHADOW  1136  12  20  0  0  0.0  12  24  0  0  0.0 

VEH  1136  HUN  NO.  TST11MAY 
SHADOW  INTERVAL  TABUIATED  BELOW 


ENTER 

LEAVE 

MO 

DY 

HR 

• 

MN 

SEC 

SYSTEM  REV 

HR 

MN 

SEC 

SYSTEM  REV 

• 

DURATION 

12 

20 

• 

• 

7 

13 

40. 

26020.2964.2 

7 

50 

32. 

• 

28232.2964.6 

2212. 

12 

20 

8 

42 

19- 

31339.2965.2 

9 

19 

11. 

33551.2965.6 

2212. 

12 

20 

10 

10 

58. 

36658.2966.2 

10 

47 

51. 

38871.2966.6 

2213. 

12 

20 

11 

39 

37- 

41977.2967.2 

12 

16 

30. 

44190.2968.6 

2213. 

12 

20 

13 

• 

8 

17. 

47297.2968.2 

13 

45 

10. 

49510.2968.6 

• 

2213. 

12 

23 

• 

• 

1 

43  16. 

6196.3009.2 

2 

20 

9- 

• 

• 

8409-3009.6 

2213. 

12 

23 

3 

11  55- 

11515.3010.2 

3 

48 

48. 

13728.3010.6 

2213. 

12 

23 

4 

40  34. 

16834.3011.2 

5 

17 

27. 

19047.3011.6 

2213. 

8.2  Test  2.  With  the  orbit  defined  in  Test  1.  the  shadow  period 
should  be  zero  around  the  autumnal  equinox.  Manual 
calculations  showed  that  the  satellite  should  be  in 
total  sunshine  from  the  first  week  in  September  to 
the  first  week  in  October. 

Function  Cards  and  Results 
*  SHADOW  1136  9  4  20  0  0.0  9  7  20  0  0.0 


VEH  1136  HUN  NO.  TST11MAY 
SHADOW  INTERVAL  TABUIATED  BELOW 


ENTER 

LEAVE 

MO  DY  HR  MN  SEC 

• 

SYSTEM  REV 

HR  MN  SEC 

SYSTEM  REV 

• 

DURATION 

9 

• 

• 

6  10  19  55. 

37195.1260.9 

10  24  33. 

• 

37473.1260.9 

278. 

9 

6  11  48  48. 

42528.1261.9 

11  52  58. 

42778.I26I.9 

250. 

9 

6  13  18  2. 

47882.1262.9 

13  21  3. 

48063.1262.9 

lfll. 

9 

6  14  47  5. 

53225.1263.9 

14  49  18. 

53358.1263.9 

133. 

REV1264.9  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 
REV1265.9  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 
REV1267.0  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 
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*  SHADOW  1136  10  6  12  0  0.0  10  8  12  0  0.0 

VEH  1136  RUN  NO.  TSTllMAY 
SHADOW  INTERVALS  TABULATED  BELOW 

ENTER  LEAVE 

MO  DY  HR  MN  SEC  SYSTEM  REV  HR  MN  SEC  SYSTEM  REV  DURATION 


REV1763.5  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 
REV1764.4  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 
REV1765.5  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 
REV1766.5  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 


10 

7  14  47  26. 

53246.1767.4 

14  51  26. 

53486.1767.5 

240. 

10 

7  16  15  6. 

58506.1768.4 

16  19  52. 

58792.1769.4 

286. 

10 

7  17  43  33. 

63813.1769.4 

17  48  43. 

64123.1769-5 

310. 

10 

7  19  ll  49. 

69109.1770.4 

19  17  44. 

69464.1770.5 

355. 

8.3  Test  3- 

An  actual  bird  observed  to  enter  sunlight 

at  a  particular 

rev.  This  orbit  was  recreated  and  the  program  was  run 

over  the  time  in  question.  This  happened  to  be  a  relatively 

eccentric  orbit.  Program  output  agreed  with  what  was  expected. 


Function  Cards  and  Results 


# 

* 

* 

* 

* 


TIME  1136  1  1  12  23  4  29  48. 

TIME  1136  1  2  12  23  4  32  48. 

TIME  1136  1  3  12  23  4  33  8. 

TIME  1136  1  4  12  23  4  38  8. 

INJKJN  1136  0  -27307.  +158.54 


RESET  TAPES  9,  10 

RESET  TAPES  2,  9 

RESET  TAPES  10,  2 

RESET  TAPES  9,  10 
+23.26  +149. 14  -922500.  +2.31  0 


*  CLR 


*  CARDS 


DEFINES  SHADOW  3000 
*  SHADOW  1136  11.  30 


lo  to  to  lo  to  So  § 


( 
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VEH  1136  RUN  NO.  TST11MAY 
SHADOW  INTERVAL  TABULATED  BELOW 


ENTER  LEAVE 


MO  DY 

HR 

MN 

SEC 

SYSTEM 

REV 

HR  MN  SEC 

SYSTEM 

REV 

DURATION 

12  24 

1 

52 

58. 

6778. 

11.4 

1  58  58. 

7138. 

11.4 

360. 

12  24 

3 

49 

13. 

13753. 

12.4 

3  54  25. 

14065. 

12.4 

312. 

12  24 

5 

45 

45. 

20745. 

13.4 

5  50  17. 

21017. 

13.4 

272. 

12  24 

7 

42 

22. 

27742. 

14.4 

7  46  6. 

27966. 

14.4 

224. 

12  24 

9 

39 

5. 

34745. 

15.4 

9  4l  1*8. 

34908. 

15.4 

163. 

12  24 

11 

36 

10. 

41770. 

16.4 

11  37  7. 

41827. 

16.4 

57- 

REV  17.4  SHADOW  PERIOD  LESS  MAN  ONE  MINUTE 
REV  18.4  SHADOW  PERIOD  LESS  THAN  ONE  MINUTE 


8.4  Test  4.  In  another  actual  case  the  time  of  entry  Into  sunlight 
was  known  (at  least  approximately)  from  telemetry. 

A  reset  tape  was  borrowed  In  order  to  test  the  program 
against  this  data.  The  times  differed  by  16  seconds. 


8.5  Test  5*  Results  from  a  similar  program  for  a  7090  computer 
were  made  available,  plus  corresponding  injection 
conditions.  The  greatest  error  in  fourteen  comparisons 
was  40  seconds.  No  systematic  bias  was  evident. 

The  reported  accuracy  of  the  7090  program  is  approximately 
one  minute. 


Function  Cards  and  Results 


*  TIME  1136  1  1  6  2  9  52  43.95 

*  TIME  1136  1  2  6  2  9  55  17.58 

*  TIME  1136  1  3  6  2  9  55  54.28 

*  TIME  1136  1  4  6  2  10  0  0.0 

*  INJRJN  1136  0  -27922.473  +8.54 


RESET  TAPES  2,  9 

RESET  TAPES  10,  2 

RESET  TAPES  9,  10 

RESET  TAPES  2,  9 
+15.421  +152.953  -1223407.5 


-.326  0 


*  CLR 


( 


*  CARDS 

DEFINES  SHADOW  3000 

*  SHADOW  1136  6  2  10  10  0.0  6  4  20  0  0.0 
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VEH  1136  RUN  NO.  TST1IMAY 
SHADOW  INTERVALS  TABULATED  BELOW 


ENTER  LEAVE 


MO  DY  HR  MN  SEC 

SYSTEM 

(7090  PCM) 

HR  MN  SEC 

SYSTEM 

(7090  POM) 

6 

2  11  52  28. 

42748. 

42720. 

12  22  28. 

44548. 

44580. 

6 

2  14  6  54. 

50814. 

50820. 

14  37  20. 

52640. 

52680. 

6 

2  16  22  . 

58920. 

58920. 

16  52  26. 

60746. 

60780. 

6 

2  IB  37  8. 

67028. 

67020. 

19  7  32. 

68852. 

68880. 

6 

2  20  52  15. 

75135. 

75120. 

21  22  38. 

76958. 

76980. 

(Note  that  7090  System  Times  are  divisible  by  60, 

,  consistent  with 

the  alleged  one-minute  program  accuracy). 


8.6  Teat  6.  In  order  to  gauge  the  effect  of  approximating  sun  position 
by  sine  curves,  a  comparison  was  made  of  two  developmental 
versions  of  the  program,  one  using  the  approximation 
technique,  the  other  using  tabulated  values  of  sun  (x,  y,  z) 
taken  directly  from  the  almanac.  Bitering  and  leaving  times 
agreed  exactly. 

Function  Cards  and  Results 
*  SHADOW  1136  6  20  20  0  0.0  6  24  0  0  0.0 

VEH  1136  RUN  NO.  TST11MAY 
SHADOW  INTERVALS  TABULATED  BELOW 


ENTER 

LEAVE 

MO  DY  HR  MN  SEC 

SYSTEM 

REV 

HR  MN  SEC 

SYSTEM 

REV 

DURATION 

• 

• 

• 

(Table  Look-up  Version) 

• 

a 

6  21  12  30  30. 

45030. 

11.7 

13  7  24. 

47244. 

12.1 

2214. 

6  21  13  59  10. 

50350. 

12.7 

14  36  4. 

52564. 

13.1 

2214. 

6  21  15  27  49. 

55669. 

13.7 

16  4  43. 

57883. 

14.1 

2214. 

6  21  16  56  29. 

60989. 

14.7 

17  33  23. 

63203. 

15.1 

2214. 

6  21  18  25  8. 

66308. 

15.7 

19  2  2. 

68522. 

16.1 

2214. 

6  21  19  53  48. 

71628. 

16.7 

20  30  42. 

73842. 

17.1 

2214. 

6  21  21  22  27. 

76947. 

17.7 

21  59  21. 

79161. 

18.I 

2214. 

‘  (Sine 

Curve  Approximation  Version) 
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m-(L) -1272/000/00 


ENTER 

MO  DY  HR  MN  SEC 

6  21  12  30  30. 
6  21  13  59  10. 
6  21  15  27  49. 

6  21  16  56  29. 
6  21  Ifl  25  8. 
6  21  19  53  48. 
6  21  21  22  27. 


SYSTEM  REV 

45030.  11.7 
50350.  12.7 

55669.  13.7 
60969.  14.7 
66308.  15.7 
71628.  16.7 
76947.  17.7 


LEAVE 

HR  MN  SBC 

SYSTEM 

13  7  24. 

47244. 

14  36  4. 

52564. 

16  4  43. 

57883. 

17  33  23. 

63203. 

19  2  2. 

68522. 

20  30  42. 

73842. 

21  59  21 

79161. 

REV 

DURATION 

12.1 

2214. 

13.1 

2214. 

14.1 

2214. 

15.1 

2214. 

16.1 

2214. 

17-1 

2214. 

18.1 

2214. 

8.7  Conclusions 

The  analysis  she  vs  favorable  results,  and  net  accuracy  at  least  as 
good  as  one  minute.  One  clue  that  It  might  be  considerably  better 
Is  afforded  by  the  DURATION,  which  during  peak  shadow  periods  does 
not  vary  by  a  second,  and  which  varies  monotanically  and  regularly 
through  those  periods  when  the  shadow  period  Is  changing  more  rapidly. 
That  is,  not  more  than  a  second  or  two  random  "Jitter"  is  noted. 
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IM-(L) -1272/000/00 


10.  0  DETAILED  FLOW  CHARTS  (Underlined  caps  are  program  tags. ) 
IO.  1  Main  Program 


EJTER 

SHADOW  Save  Index  Registers 

Read  ^  file  from  Reset  Tape 

|r 

Print  Vehicle  and  Run  numbers,  and  headers 

Rev  no.  input?  -  Yes - WJsing  LTT,  Set  WHISTLE  and  Start  times  from 

I  rev  number 

NO 

Set  Universal  time,  get  machine  time 

Jr 

Set  WHISTLE  »  quitting  time 
Set  TAU  »  start  time 


I  X - 

HKEEP  Clear  TAUTOTAL,  Compute  PER/ 8 
v 

Set  TAUSTAHT  -  TAU 

Compute  XBASE,  YBASE  for  this  year 

SEARCH  Set  XR  6  -  0 
-  1 

In  sunshine  this  t?  -  NO-  ■■  »Get  XR  6  ■  1  for  two  passes  through  SEARCH 
1  to  get  a  better  TAUIN 

YES 

i 


INCR3  Increase  TAU  by  one  minute  to  establish 

direction  (toward  or  away  from  x',  y',  z') 

In  sunshine  this  t?  -  NO - 


to 


STCHAD 


yis 


to 


INCR 
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incr 


Increase  TAU  by  l/8  period 

Check  for  < one  (DONECHCK  Subroutine) 

y 

In  sunshine  this  t?  -  NO - 

YES 

Are  we  getting  closer  to  x1,  y1,  z'  ? 
YES  1 - NO 


r 

Had  we  been  getting  closer? 

- io  YES  r- 

FINER  Decrease  TAU  by  PER/8 
Set  DIFF2  yla  COMPUTE 


- 1  r 

INCR1  Store Tmodiried  TAU 


tills 


t?  -  N0- 


Oetting  closer?  -  NO - WRNODATA 

(Subroutine) 


INCR4  Increase  TAU  by 
One  Minute 


to  SEARCH 


f 

ate 


STSHAD  Save  entering  times 

INCR2  1  Increase  TAU  by  one  minute 

(First  time,  step  through  shade 
in  one  minute  increments) 

l 

Sunshine  this  t? 

_ ,  pF3 

Set  TAUOUT  -  TAU 

1 

Print  one  line  of  data 

Did  we  start,  in  shadow? 

r— YES  r-NO 

to  INCR  to  CYCLE 
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1*4- (L)  -1272/000/00 


CYCLE  Set  TAU  -  TAUIN  +  period  -  one  minute 
Check  for  done  (DONECHECK) 

i  . 

Sunshine  this  t?  -  NO - 

I 

YES 

1 


INCRE 


Increase  TAU  by  one  minute 
Sunshine  this  t?  -  NO - 


Y^S 


Are  we  getting  closer  to  (x'.y ' ,z' )1 

i  I 

YES  NO 

_ 1  1 

PRNODATA  Subroutine 

i 

to  SEARCH 


SHAKE 


Decrease  by  one  minute 

i 

Sunshine  this  t? 

i  A 

NO  YES 

I 


INCRElr->Increase  TAU  by  one  second 
Sunshine?  -  NO 

i  I 

YES  Save  entering  times 
-(SAVE  subroutine) 


• - 1  r - suDroutme j 

Set  TAU  -  TAUIN  pluS  last  DURATION  (TAUDUR) 


INCRE2  Store  modified  TAU 

“-i 


In  sunshine  this  t?  -  NC 

YES  SHADE1  (—►Increase  TAU  by  one  second 

l  I  I 

In  sunshine  this  t? 

I  i 

NO  YES 

l 


Decrease  TAU  by 
one  second 


Print  one  line  of  data 

I 

to  CYCLE 
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DETAILED  FLOW  CHARTS  (Cont'd) 

10.2  COMPUTE  SUBROUTINE 

The  inputs  to  this  subroutine  are  TAU,  X,  Y,  Z,  and  the  inertial  coordinates  of 
the  sun.  The  outputs  are: 

p 

UCDSQ  *  D,  where  D  is  the  distance  from  (x',y',z')  to  the  satellite. 

2 

LCDSQ  -  d  ,  where  d  is  the  distance  from  (x',y' ,z*)  to  the  border  of  the 
shadow. 

DIFFNOW  *  (D^  -  d2).  If  negative,  satellite  is  in  shadow. 

DIFFLAST  -  (D  -  d  )  for  the  previous  entry  to  this  subroutine. 

p  p 

DIFF2  *  (DIFFNOW  -  DIFFLAST),  or  (D  -  D  -)  for  two  successive  subroutine 
executions.  n  n" 

LSTDIFF2  =  Previous  value  of  DIFF2. 

(These  "second  differences"  are  used  in  the  SEARCH  portion  of  the  program,  which 
locates  initial  shadow  boundaries.  They  are  also  essential  to  the  logic  which 
determines  that,  for  a  given  rev,  shadow  interval  less  than  60  seconds.  ) 

Enter  Subroutine 

i 

Commute  Universal  Time 
Set  DIFFLAST  and  LSTDIFF2 

l  ,  x 

Commute  elements  this  TAU  (TTE) 

Convert  to  rectangular  (PTR) 

Confute  RHOSQ,  radius  squared  of  sphere  containing  satellite. 

Compute  time  of  year  in  hours,  call  it  HANGLE 
HANGLE  <  XBASE?  -  YES  * AA10  Divide  by  XPERN 

i 


NO 
1 

HANGLE  >  XBASE  +  XPER  ?  YES- 
NO 

Divide  bv  XPER _ 


I 


c 


BB10  Divide  by  XPERN1 


1 


Multiply  by  IT  to 
convert  to  radians 

I  9 - 

Store  as  XANGLE 
for  checkout  purposes 

i 


Multiply  by  II 

Add  II,  converting  HAJTCLE 
to,  radians 
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I 

Fetch  Cosine  this  angle 

1 

Subtract  fudge  constant 

i 

Convert  to  feet 

i 

Store  cm  SUNX 

[ 

Accumulate  range  of  sun  (squared) 

1 

Similarly,  from  HANGIE  compute  YZANGLE  in  radians 

1 

Fetch  sine  this  angle 

i 

Negative  ?  -  YES  Normalize  by  YNORMN 

io 

i 

Normalize  by  YNOHMP 


i  { - 

Convert  to  feet 

Store  as  SUNY 

1 

Accumulate  range  of  Sun  (squared) 

Similarly,  compute  SUNZ 
Commute  range  of  sun 
Commute  x',y',z' 

Commute  UCDSQ  and  LCDSQ 

Set  DIFFNOW,  DIFF2,  leave  DIFFNOW  in  accumulator 


Exi 


i 
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DETAILED  FLOW  CHARTS  (cont'd) 


10. 3  PRDATA  Subroutine 

This  subroutine  Is  entered  whenever  the  satellite  leaves  shadow.  It  prints  on- 
and  off-line  the  times  of  entering  and  leaving,  the  latitude/ longitude  of  the 
subpoint  at  these  times,  and  the  duration  of  time  in  shadow. 


Enter  Subroutine 

4 

Set  up  output 

4 

Print  on-line 

4 

Print  off-line 
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DETAILED  FLOW  CHARTS  (cont'd) 

10.4  SAVE  Subroutine 

This  subroutine  saves  time  and  position  vhen  the  satellite  enters  shadow. 

Enter  Subroutine 

Store  desired  values  of  time 

I 

Set  up  for  LATLONG  subroutine 

Compute  and  store  latitude  and  longitude  for  printout 
Exit 
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DETAILED  FLOW  CHARTS  (cont'd) 


10. 5  DONECHCK  Subroutine 


This  subroutine  Is  entered  when  IAU  Is  incremented,  to  check  for  "done", 
so,  a  Summary  printout  Is  presented  showing  per  cent  of  time  In  shadow. 


Enter  Subroutine 

t 

TAU  <  WHISTLE?  -  YES - >Exit 

NO 

i 

Print  per  cent  time  in  shadow 

1 

Reset  Index  Registers 

1 

Return  to  COP 


If 
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DETAILED  FLOW  CHARTS  (Con't.) 

10.6  DECAY  SUBROUTINE 

This  subroutine  Is  entered  vhen  RAD ^  Earth  radius,  Implying  that  the  bird 
has  come  down. 

Enter  Subroutine 

Print  on-line  ORBIT  HAS  DECAYED 

i 

Print  summary,  reset  Index  registers 
Return  to  COP 
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DETAILED  FLOW  CHARTS  (Coat'd) 

10.7  IATLONG  SUBROUTINE 

This  subroutine  computes  geodetic  latitude  suid  longitude  in  degrees  and  minutes 
and  sets  up  one  part  of  a  two-part  table,  corresponding  to  entering  and  leaving, 
depending  upon  index  register  U,  which  is  preset  to  zero  for  entering  and  six  for 
leaving. 


filter  Subroutine 


T 


Set  Character  for  North  or  South  latitude 
Cot  vert  to  geodetic,  degrees  and  minutes 
Stcre  in  table  for  output 

1  character  for  East  or  Vest  longitude 
ert  to  degrees,  minutes 
Store  in  table  for  output 


22  May  1963 


32 


tm-(l) -1272/000/00 


DETAILED  FLOW  CHARTS  ( Coat'd) 

10.8  FRNOSHAD  SUBROUTINE 

This  subroutine  Is  entered  whenever  the  shadow  period  was  less  than  sixty 
seconds  for  a  given  rev.  It  prints  an  appropriate  comment  on-line. 


filter 

i 


Print 


Subroutine 

REV  XX.  SHADOW  PERIOD  LESS  THAN  CHE  MINUTE 


I 
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Reports  that  SHADOW  Is  a  1604  COP- 
cont rolled  program  to  compute  the 
tlmea  at  which  a  satellite  enters 
and  leave a  the  shadow  of  the  earth. 
Alao  reports  that  results  are 
printed  on-line  and  off-line. 
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